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. . . The interpretation of results of recent r decay determinations of | Vus \ , which yield 

^ \ values ~ 3(T low compared to 3-family unitarity expectations, is complicated by the 

\ slow convergence of the relevant integrated D = 2 OPE series. We introduce a class 

' of new sum rules involving both electroproduction and r decay data designed to 

^ ■ deal with this problem by strongly suppressing D = 2 OPE contributions at the 

correlator level. Experimental complications are briefly discussed and an example 
! of the improved control over theoretical errors presented. The uncertainty on the 

resulting determination, \Vus\ = 0.2202(39), is entirely dominated by experimental 
1-^^' errors, and should be subject to significant near-term improvement. 

JiLi' PACS numbers: 12.15.Hh,13.35.Dx,11.55.Hx 
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■ I. INTRODUCTION 

(N 

> , , 

O ! The CKM matrix element, |Kts|, is one of the fundamental parameters of the Standard 

Model (SM). Determinations from multiple sources can help to improve the accuracy with 
which it is known and/or test for the presence of non-SM contributions in strangeness- 
changing weak processes. Current analyses of Ki^ and r[i^^2]/r[7r^2] [H, using lattice 
input for /+(0) and fx/ f-K, respectively jl, S], yield values which are in good mutual 
agreement and compatible with the expectations of 3-family unitarity jl, 0, 0]- lYyA 
can also be obtained using flavor-breaking (FB) hadronic-r-decay-based sum rules SUS, 
lol . ITl| . Recent r decay analyses [l^, 11], employing updated information 12. 13. uA 



in 



oo 
O 



X 



isl . [igI . [itI . 18| on the older strange decay distribution [H, 2^ , yield values ~ 3a below 
c3 ' 3-family-unitarity expections. 

In existing r decay determinations, the dominant source of error on \Vus\ is the un- 
certainties on weighted integrals over the inclusive strange decay distribution. This error 
will be signiflcantly reduced by ongoing B-factory analyses. Nominal theoretical errors, 
estimated with conventional prescriptions, are small, holding out the prospect of results 
competitive with those from K^^ and r[i^'^2]/r[7r^2]; once the B-factory analyses are 
complete. A potential complication, however, arises from the slow convergence of the 
relevant integrated D = 2 OPE series. Evidence suggests that theoretical errors may be 
underestimated (in some cases, signiflcantly) as a consequence of this behavior. 

In this paper we consider alternate sum rules for |V^s|, involving both r and electropro- 
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duction, rather than just r, spectral data. The combinations chosen have, by construc- 
tion, already at the correlator level, a strong suppression of the potentially problematic 
D = 2 OPE series, and hence also a strongly reduced D = 2 trunctation contribution to 
the theoretical uncertainty. The rest of the paper is organized as follows. In Section [ITl we 
first briefly outline the purely r decay approach and associated D = 2 OPE convergence 
problem. Then, in Section IIIII we introduce and discuss the alternate, mixed r decay- 
electroproduction sum rules. Finally, Section HVl outlines the spectral and OPE input, 
discusses some experimental complications, and provides an illustration of the utility of 
the mixed sum rule approach. 



II. THE HADRONIC r DECAY DETERMINATION OF 

For any correlator, 11, without kinematic singularities, and any analytic weight, w{s), 
analyticity implies the finite energy sum rule (FESR) relation, 

so I r 

w{s) p{s) ds = (b w{s) n(s) ds , (1) 

27rz 

where p{s) is the spectral function of n(s) and the OPE expansion of n(s) can be 
employed on the RHS for sufficiently large Sq. \Vus\ is obtained by applying this relation 

to the FB correlator difference AIlT-fs) = 



n(o+i) (s) - n^°+^^ (s) 



where Hy/^.j^ are 



the spin J = 0, 1 components of the flavor ij, vector (V) or axial vector (A) current 
two-point functions, and the corresponding spectral functions, Pyjy^.ij, are related to the 
differential distributions, dRv/A;ij/ds, of the normalized flavor ij V or A current induced 
decay widths, Rv/A;ij = r[r~ — > 1/^- hadionsy /A;ij (7)]/r[T~ 1/^-6^ z/e (7)], by j2l| 



dRv/A;ij _ 12TT'^\Vij\'^SEW 



ds m 
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V/A;ij\'^> "^L \yTJh'V/A;ij^ 



with = s/ml, w^^'°\y) = {l-y)\l + 2y), w^^'°\y) = 2y{l-yY, Vij the flavor CKM 
matrix element. Sew a short- distance electroweak correction |22|, and (0 -|- 1) denoting 
the sum of J = and 1 contributions. The J = contribution to any J = + 1/J = 
decomposition will be referred to as "longitudinal" in what follows. 

The use of the J = -|- 1 difference. An,-, is the result of the extremely bad behavior 
of the integrated longitudinal D = 2 OPE series [23], which precludes working with FB 
FESRs based on the linear combination of J = 0, -|- 1 spectral functions appearing in 
Eq. ([2]). The Pv^A-iji^)-: and from these, Ap^ are obtained by identifying and subtracting, 
bin-by-bin, the longitudinal contributions to dRy/A-ijI ds. This can be done with good 

accuracy because, apart from the vr contribution to p^A-ud ^"^^ ^ contribution to p^a-^s 
(which are determined by /^r and /x, respectively, and hence very accurately known) all 
contributions to PyjA-ij doubly chirally suppressed, by factors of 0[(mj =F '^j)^]- The 
ij = ud longitudinal contributions are thus, to high accuracy, saturated by the tt pole 
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term. Continuum longitudinal ij = us contributions, which are small, but not entirely 
negligible, are determined from dispersive [1^ and sum rule 2^ analyses of the strange 
scalar and pseudoscalar channels, respectively, analyses which are strongly constrained 
by their implications for rus 2^ . 

Given w{s) and sq < m^, \Vus\ is determined by first constructing, from the longitu- 
dinally subtracted dRv/A;ij/ds, the spectral integrals 



so 



ds 



mt 



and, from these, the FB combinations. 



(3) 



(4) 



Using the OPE representation of 5i?y_,_^(so), and inputting \Vud\ and the required OPE 
parameters from other sources, one obtains |8|, from Eq. ([T]), 



IK, 



R 



V+A;us 



So) 



R 



V + A;u 



- [SR^_,A^o)] 



OPE 



(5) 



Since 
R 



at scales Sq 



3 GeV' 



6R 



V/A 



So) 



J OPE 



is typically much smaller than 



V/A;ud,us 



[sq) (usually at the few-to-several- % level), Eq. ([5]) yields a determination of 



\Vus\ with a fractional uncertainty much smaller than that on 



SR^/Ai^o) 



itself M. 



OPE 



A particularly advantageous case, from the point of view of experimental errors, is 



that based on sq 



m: 



and the weight w{s) 



W^(00)(l/r) = Wj, 



(0,0) 



(yr) In this case, the us 



and ud spectral integrals appearing in Eq. ([5]) are fixed by the total strange and non- 
strange r branching fractions, allowing one to take advantage of improvements in the 
errors on a number of the strange branching fractions in advance of the completion of the 
remeasurement of the full us spectral distribution. A disadvantage of this approach is 
that, working with only a single Sq, one is unable to test the stability of the output |\4s| 
values with respect to sq, a crucial step to ensuring that estimates of the accompanying 
theoretical uncertainty (which, in some places in the literature, are quoted to be as low 
as 0.0005) are sufficiently conservative. See below for more on this point. 

The OPE representation of 6Ry^j^{so) is, of necessity, truncated, in both dimension 
and the perturbative order of the relevant Wilson coefficients. Estimating the associated 
theoretical uncertainty is complicated by the less-than-ideal convergence of the J = + 1, 



D = 2 OPE series. Exphcitly [27|, [28 



[An.(g 



2^-]OPE 
D=2 



3 ■ms{Q'^ 



27r2 



1 + + 19.93a2 + 208.750^ + 



(6) 



with a = as{Q^)/n, and as{Q^) and ms{Q^) the running coupling and strange quark 
mass in the MS scheme. Since independent determinations of as [2I, [13, 31, HI, [H, 



4 



34 . [ssl . 36| imply a(m^) ~ 0.10, the convergence of the D = 2, J = + 1 series at 
the spacehke point on the OPE contour is marginal at best. While (at least if one 
works with the contour improved (CIPT) prescription 37|, in which the large logs are 



resummed point-by-point along the contour) the convergence of the integrated series can 
be improved through appropriate weight choices [38|, taking into account that |a<j(Q^)| 
decreases as one moves away from the spacelike point along the contour, one expects, in 
general, rather slow convergence, which makes conventional truncation error estimates 
potentially unreliable. Fortunately, the growth of with decreasing sq means that 
omitted higher-order terms become relatively more important at lower scales, and hence 
that any premature truncation of the slowly converging integrated D = 2, J = + 1 series 
will show up as an unphysical so-dependence in the extracted values of | Vus \ ■ Unphysical 
So-dependence can also be produced by incorrect input for poorly known, or unknown, 
condensates relevant to D > 4 OPE contributions (-0 = 6 and 8 in the case of W(oo))- 

Such unphysical so-dependence is, in fact, seen, at a scale significantly larger than the 
estimated D = 2 truncation error, in recent r decay analyses [llj. This is illustrated in 
Figure [H which shows results for W(oo), and three additional weights, Wiq, Wiq and W20, 
introduced originally to improve the integrated D = 2, J = + 1 convergence [38i]. Of 
particular note is the situation for the experimentally favorable W(oo) weight case, where 
the instability is much larger than full estimated theoretical error. 

It is worth noting that the D = 2 truncation component of the 0.0005 theoretical error 
in the W(oo) case is obtained by combining an uncertainty associated with the residual 
scale dependence with the shift obtained by omitting the last term included in the D = 2 
series, all evaluations being performed using the CIPT prescription and the truncated 
D = 2 Adler function form. Alternate evalutions, using the truncated correlator (rather 
than truncated Adler function) form, and/or using the fixed order (FOPT) rather than 
CIPT prescription, are, however, also possible. At a given, common truncation order, 
all such evaluations are equivalent to the CIPT Adler function evaluation, differing from 
it only by corrections of yet-higher order. While the difference of \Vus\ values obtained 
using the O(a^) and O(a^) CIPT Adler function evaluations is, indeed, small (5|V^us| = 
—0.0003), shifting to alternate D = 2 evaluation schemes leads to much larger shifts. 
For example, shifting from the 0{a^) CIPT Adler function evaluation to the 0(a'^) CIPT 
correlator version yields instead 5|V^us| = —0.0008, while shifting from the O(a^) CIPT 
Adler function (correlator) versions to the 0(a^) FOPT correlator version yields the 
even larger shifts 6\Vus\ = 0.0019(0.0023) |39|]. With plausible arguments in favor of 



both the CIPT and FOPT prescriptions in the literature [31|, l37[, such shifts suggest 
the conventional D = 2 truncation error estimate, which leads to the total estimated 
theoretical uncertainty, 6\Vus\ = 0.0005, for the sq = m^, W(oo) determination, is far from 
a conservative one. 

In view of the possibility of much-larger-than-anticipated D = 2 truncation uncer- 
tainties on the values of \Vus\ extracted using the An,- FESRs, we consider, in what 
follows, FESRs based on alternate correlator differences designed to have, already at the 
correlator level, much reduced D = 2 contributions. Such FESRs also allow one to in- 
vestigate whether the sizeable So-instability observed in the results of the iy(oo)-weighted 
AIIt- analysis is a consequence oi D = 2 truncation uncertainties, or of unexpectedly 
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large D = 6,8 OPE contributions. 



III. NEW MIXED r-ELECTROPRODUCTION SUMRULES FOR |K 



Problems associated with the slow convergence of the integrated D = 2, J = 0+1 OPE 
series can be reduced by considering alternate FESRs based on correlator differences, 
An, sharing with AII^ the vanishing of D = OPE contributions, but having D = 2 
contributions suppressed at the correlator level. Since a V/A separation of the flavor us 
decay distribution is not presently feasible. All should involve the us V+A combination. 
The leading order term in the D = 2 Wilson coefficient can be removed by forming the 
appropriate difference of Tfy+A-us electromagnetic (EM) correlator. Hem- The 

following combinations (having the same Hy^^^.^g contribution as AII^) have, in addition, 
vanishing D = contributions: 



An, ^ gn^M - n^vtXus - 2(2 + ^)u)^;j + 2^u'X 



(0+1) 



r(o+i) 



ud 



(7) 



The K = 1/2 combination is strictly FB. Bearing in mind that a(m^) ~ 0.1, the corre- 
sponding D = 2 OPE contribution. 



[An«(g' 



OPE 3 ms{Q 



2\ r 



^=2 27r2 g2 



-a + 4.3839a^ + 44.943a^ + 



is seen to be strongly suppressed, by more than an order of magnitude, compared to 
[An^-j^p^. A similar suppression turns out to be operative for the D = 4 contributions. 
Explicitly, up to numerically tiny O(m^) corrections, and neg lecting, for simplicity of 
presentation, r = (m^ — mu)/{ma + my), one has, to 0{d?) jU 4o| . 



[An.(g2)] 
[An«(g')] 



2^1 OPE 
Z)=4 



2^1 OPE 
D=4 



2 

2 



13 

[{miii) — (rrisSs)] I 1 — a a" 



4- 16k 



a + 



59 - 236k 
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+ -a H a (msSs) 

'3 6 ' \ ^ / 



(9) 



(10) 



where in both cases the strange condensate term is numerically dominant. 

Defining RIm{so) = [1271^ Sew/ ml] j^' dsw{s)pEM{s) and [5R:{so)f'''' = 
\i2'K'^ S Ew / m^] §^^^^^^^dsw{s) [/\Tlt^{syf^^ , one then has, for any analytic w{s) and any 
So large enough the OPE representation is reliable, the An, analogue of Eq. ([5]), 



IK. 



^V+A-usi-^o) 



\ QReaAso) 



2(2+K)i?™.„^(so)-2«;i?« 



OPE 



fill 
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The suppression of the D = 2 and 4 contributions in [An^]'^'^^ does not persist to 
higher D. For example, with rc = (ss) / {ii), the D = 6 contributions, in the vacuum 
saturation approximation (VSA), become 



OPE _ vra 

D=6;VSA ~ Q6 



81 ^ 



\AU (0^)1 = 

[l^l-l-nW )\d=6;VSA Q6 



-32 + 128k\ 32r 



.21 



9/9 



(12) 
(13) 



m 



typically significantly larger for AIIk than for AII^. 

To deal with such potentially enhanced, but phenomenologically poorly determined 
D > 4 contributions, it is useful to employ polynomial weights, w{y) = ^^=0^™?/ 
with y = s/sq. Integrated D = 2k + 2 OPE contributions then scale as 1/sq. The 
strong suppression of D = 2,4 contributions, which scale more slowly with sq, then 
means one can, for example, employ the VSA estimate for D = 6, and ignore D > 6 
contributions, but look for sq values large enough that \Vus\ becomes stable with respect 
to So, indicating that D > 4 contributions and/or deviations from the input assumptions 
about these contributions have decreased to a negligible level. 

The expected enhanced role of D = 6 and higher contributions in An^ means that 
higher degree weights like wio, W20 and wio, introduced to improve the integrated D = 
2 convergence for the An,- FESRs, are likely to represent less useful choices for the 
AIIk analysis. The strong suppression of D = 2 contributions, however, opens up the 
possibility of using weights which provide less good integrated D = 2 convergence but 
better control over integrated D > 4 contributions. Thus, e.g., if it is the slow D = 2 
convergence which is responsible for the significant so-instability of the W(oo)-weighted 
AIIt- FESR results shown in Figured! the analogous An^ FESR might be rendered stable 
by the reduced D = 2 contributions, allowing improvements in the strange branching 
fractions (whose sum provides an improved determination of Ry^^A-usi^D) ^e used in 
reducing the error on the numerator in Eq. (fTT]) for w = W(oo) and Sq = Tn"^. Similarly, 

it might become possible to employ the weights, WN{y) = 1 — -j^zjU + if^' which, like 
1(7(00), display slow integrated D = 2 convergence for AIIt- but are useful for handling 
D > 4 contributions (written generically as X]d=6 8 ■•• ^-d/Q^) since (up to corrections 
of 0{[as{ml)f)) only a single integrated D > 4 contribution, [—1)'^C2n+2/[{N — l)s^], 
survives on the OPE side of the wn FESR. 



IV. INPUT, COMPLICATIONS, RESULTS, AND DISCUSSION 

In this section we illustrate the utility of the new mixed FESRs and point out some 
experimental complications, focussing on the An«;=i/2 case, whose D > 4 contributions 
vanish in the SU{3)f limit, and are thus expected to be optimally suppressed. 
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FIG. 1: \Vus\ versus sq from the AH^ FESRs 
for t(72o, wiQ, wio and W(oo)- 



FIG. 2: \Vus\ versus sq from the An^=;^/2 
FESRs for t«(oo)> ^^^lOj and w^. 



A. OPE input 

To suppress OPE-breaking contributions from the region of the contour on the RHS 
of Eq. (pP) near the timehke point on the contour, we restrict our attention to w{s) having 
a zero of order > 2 at s = Sq and to Sq > 2 GeV^ |i41i] . 

D = 2 OPE integrals are evaluated using Eq. and the CIPT prescription j37| , 
with as(Q^) and ms{Q^) the exact solutions associated with the 4-loop-truncated /3 and 
7 functions [i^ and the initial conditions, ms(m^) = 100 ± 10 MeV fio'], as(m^) = 
0.323(17). The latter is obtained by running a very conservative assessment, 0.1190(20'), 
of the average of several recent independent as(M|) determinations [29 



33, m 



36j down to the r scale using the standard self-consistent combination of 
4- loop running and 3-loop matching at the flavor thresholds 43|. To be conservative, 
we assign the sum of absolute values of the contributions of all computed orders as 
the truncation component of the D = 2 uncertainty (producing a 100% uncertainty 
if all contributions have the same sign, larger otherwise). The error on the truncated 
D = 2 sum associated with that on the overall [ms(m^)]^ factor is also evaluated using 
the conservative all-absolute-values prescription. The truncation and m^-scale errors 
are combined in quadrature with the much smaller error induced by the uncertainty on 
as(m^) to obtain the full D = 2 error. 

D = A OPE input and uncertainties are as follows, {meii) is fixed using the 
GMOR relation, (rrisss) using conventional ChPT quark mass ratios 4J] and the value, 
Tc = {ss) / {ii) = 1.2 ± 0.3 obtaiiied by updating the quenched-lattice-data-based deter- 
mination, Tc = 0.8(3), of Ref. |45j] using the average, fss/ fs = 1.21(4) [46], of recent 
Uf = 2 + 1 lattice determinations j47| . The strange condensate term dominates both the 
D = 4 contribution and its error, but produces only a very small impact on \Vus\ as a 
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consequnce of the suppression of the coefficient function seen in Eq. (fTU]) . 

D > 4 contributions involve poorly known or phenomenologically undetermined con- 
densate combinations. We estimate D = 6 contributions using the VSA and ignore D > 8 
contributions. If D > 4 contributions are small, the details of these assumptions are ir- 
relevant. If not, and the assumptions are inaccurate, the l/sg (1/sq, ■ ■ ■) dependence of 
integrated D = 6 (8, ■ ■ ■) contributions will lead to an unphysical so-dependence of \Vus\- 
We look for weights which produce a good window of So-stability in order to ensure that 
D > A contributions are either negligible or estimated with sufficient accuracy. 



B. Spectral input 



Results for Ry/A-udi^o) Rv+A;usi^o) are based on the ALEPH us [l9^ and ud 48 
spectral data, for which information on the relevant covariance matrices is publicly avail- 
able. The ud data has been modified to incorporate the recent improved V/A separation 



for the KKn mode [30|] made possible by the BaBar determination of the 7 = 1 KKn 
electroproduction cross- sect ions [i^. A small rescaling is applied to the continuum ud 
V+A distribution to reflect changes in Sew, Rv+A;us and the electron branching frac- 
tion. Be- With the lepton-universality-constrained result = 0.17818(32) [sT] and an 
updated total strange branching fraction Bus = 0.02858(71), the ud normalization is 
Rv+A;ud = 3.478(11). For \Vud\, we use the latest update, 0.97425(23), from the super- 
allowed nuclear 0+ ^ O"*" /? decay analysis 0]. 

Though BaBar and Belle have not completed their re-measurements of the inclusive us 
distribution, dRv+A;us/ds, an interim partial update can be obtained (following Ref. |5Q] ) 
by rescaling the 1999 ALEPH distribution [19], mode-by-mode, by the ratio of new to old 



branching fractions. The new branching fraction results are taken from Refs. [12, ll3|, [1 



15|, Il6|, [171, Il8| . Unfortunately, this strategy does not allow the corresponding covariance 
matrix to be updated. The improved precision on the new strange branching fractions 
can thus be translated into a correspondingly improved us spectral integral error only for 
w = W(oo) and sq = m^. Since the recently measured K and vr branching fractions [ll] 
are compatible with SM expectations at the ~ 2a level, we evaluate the vr and K pole 
spectral integral contributions using the more precisely determined 7r^2 and Kf^2 input. 

REAii^o) is obtained from the EM spectral function, Pem{s), which is related to the 
bare e+e~ —>■ hadrons cross-sections, ai,are{.s), by 

Pem{s) = s abare{s) /l67r^aEM{0Y . (14) 

It is well known that problems exist with the compatibility of the measured vrvr and 
7r+7r^7r'^7r° cross-sections and those implied by / = 1 r decay data, even after known 
isospin-breaking corrections are taken into account [52|. Preliminary BaBar n^TT'TT^ir^ 
cross-section results [53f| reduce considerably the latter discrepancy, but have not yet 
been flnalized. The situation for vrvr is somewhat muddier. The recent KLOE update [5J] 
yields results now in reasonable agreement with CMD2 and SND below the p peak and 
with a reduced discrepancy above it, while preliminary BaBar results [55|] are instead 
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in better agreement with r expectations. In addition, the most recent r-based analy- 
sis j32| produces an as{Mz) in excellent agreement with two recent high-precision lattice 
determinations [33|, while electroproduction-based analyses (albeit not updated for new 
post-2005 experimental results, and without the careful fitting of D > 4 OPE contri- 
butions performed for the r case) yield values ~ 2a too low |56|, again favoring the r 
version of the 1=1 spectral distribution. We thus deal with the 1 = 1 discrepancies 
by replacing both tttt and Att EM results with the corresponding r expectations. Since 
(i) the V/A separation for the KKir contribution to r decay has been performed using 
CVC and the BaBar 1 = 1 EM cross-sections and (ii) the tttt, Att and KKn contributions 
largely saturate Pv^^ii^) below s = ml, this is effectively equivalent to replacing Ani/2 
with the alternate combination 

where Ilv-i=o is the 1 = octet analogue of ^y^d ■ cross-sections are taken from 
Whalley's 2003 compilation js^] and recent updates reported in Refs. 58, 59, 6^. Where 
needed, vacuum polarization corrections are computed using F. Jegerlehner's code 61 



C. Results and discussion 

The results for \Vus\ as a function of sq obtained from the Ani/2 FESRs for 1^(00)5 ^2, 
W3, W4 and the weight, wio, producing the best AH,- so-stability plateau in Figured] are 
displayed in Figure [2l In all but the last case a very good so-stability plateau is found. In 
addition, the \Vus\ obtained at the highest accessible scale, Sq = (the right endpoints 
of the curves) are all, without exception, in extremely good agreement. The very good 
stability plateau for if(oo) strongly suggests that the instability seen in the analogous 
AIIt- analysis was a result of the slow D = 2 convergence. In contrast, the quality of the 
stabihty plateau for wio has deteriorated in going from An,- to Ani/2, most likely due 
to the increased size of D > 4 contributions. Even so, the \Vus\ values for wio converge 
nicely to the stable results from the other weight cases as sq ml. 

Because of the very good stability found for W(oo), it is possible to quote a final 
determination based on w = W(oo) and sq = m^, a choice which allows us to benefit from 
the improved BaBar and Belle strange branching fraction determinations. We find 

iKsl = 0.2202(27)„,(28)£;Af(2)y;„rf(4)^.„rf(2)oPs = 0.2202(39) (16) 

where the errors are those associated with the inclusive us branching fraction, the residual 
1 = EM spectral integral, the residual inclusive ud V and ud A branching fractions, 
and the combined D = 2 and D = 4 OPE contribution, respectively. 

While, within current errors, the result for \Vus\ is compatible with either 3-family- 
unitarity or the recent AH,- determinations, and hence does not help in resolving the 
~ 3a discrepancy between the two, prospects exist for siginificantly reducing the main 
components of the error. First, errors on the weighted / = EM integrals will be reduced 
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through ongoing work on the exclusive EM cross-sections at VEPP2000, BaBar and 
Belle. Second, errors on the us spectral integrals will be significantly reduced by BaBar 
and Belle analyses of both the branching fractions of as-yet-unremeasured strange modes 
(including the sizable K^ir^n' and previously estimated, but unmeasured, K3it and KAtt 
modes) and the inclusive us V+A distribution. Obtaining the inclusive us distribution, 
and not just the branching fractions, is crucial to performing the So-stability checks, 
themselves crucial to demonstrating that D = 2 convergence and D > 4 contributions 
have, indeed, been brought under good control. To reduce the us-distribution-induced 
contribution to the error on \ Vus\ to, e.g., the ~ 0.0005 level requires ~ 1.3 x 10~^ precision 
on the inclusive us branching fraction, and hence, almost certainly, pursuing previously 
undetected higher multiplicity modes having branching fractions down to the few xlO~^ 
level. 

We close by stressing the complementarity of the AII^ and Ani/2 analyses. The latter, 
by construction, has significantly reduced OPE-induced uncertainties. The smallness of 
the OPE contributions to the denominator of Eq. ([5]), however, means that global normal- 
ization uncertainties common to the ud and us spectral distributions cancel, essentially 
entirely, in the AII^ determination. This is not the case for the Ani/2 analysis, where EM 
and T normalization uncertainties are independent, leading to an increased experimental 
error on |\4s|- As we have seen already in the W(oo) case, employing the same weight in 
both FESRs and comparing the so-dependences of the resulting \ Vus\ determinations can 
also help in shedding light on the source of any so-instabilities found in the AII^ analysis, 
where OPE-induced errors are more difficult to reliably quantify. 
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